Mon. Not. R. Astron. Soc. 000, ??-|8](2011) 



Printed 19 January 2013 



(MN I^^T^ style file v2.2) 



The haloes of bright satellite galaxies in a warm dark 
matter universe 



O 



Mark R. Lovell^*, Vincent Eke\ Carlos S. Frenk^, Liang Gao^'^, Adrian Jenkins^ 
Tom Theuns^'^, Jie Wang^, Simon D. M. White*^, Alexey Boyarsky^'^, 
and Oleg Ruchayskiy 

^ Institute for Computational Cosmology, Durham University, South Road, Durham, UK, DHl 3LE 
'^National Astronomical Observatories, Chinese Academy of Science, Beijing, 100012, China 
'"^Department of Physics, University of Antwerp, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium 

Max-Planck-Institut fiir Astrophysik, Karl-Schwarzschild-Strafie 1, 85740 Garching bei Miinchen, Cermany 
^ Ecole Polytechnique Federale de Lausanne, FSB/ITP/LPPC, BSP 720, CH- 101 5, Lausanne, Switzerland 
^ Bogolyubov Institute of Theoretical Physics, Kiev 03680, Ukraine 

CERN Physics Department, Theory Division, CII-1211 Ceneva 23, Switzerland 



o 
u 

o 

(N 
> 

(N 

a^ 
o 



Accepted 2011 November 13. Received 2011 October 21; in original form 2011 April 15 



ABSTRACT 

High resolution N-body simulations of galactic cold dark matter haloes indicate that 
we should expect to find a few satellite galaxies around the Milky Way whose haloes 
have a maximum circular velocity in excess of 40 kms^^. Yet, with the exception 
of the Magellanic Clouds and the Sagittarius dwarf, which likely reside in subhaloes 
with significantly larger velocities than this, the bright satellites of the Milky Way 
all appear to reside in subhaloes with maximum circular velocities below 40 kms^^. 
As recently highlighted by Boylan-Kolchin et al., this discrepancy implies that the 
majority of the most massive subhaloes within a cold dark matter galactic halo are 
too concentrated to be consistent with the kinematic data for the bright Milky Way 
satellites. Here we show that no such discrepancy exists if haloes are made of warm, 
rather than cold dark matter because these haloes are less concentrated on account 
of their typically later formation epochs. Warm dark matter is one of several possible 
explanations for the observed kinematics of the satellites. 
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1 INTRODUCTION 

Measurements of temperat ure anisotropies in the microwave 
background radiation (e.g. iKomatsu et al. 2011 ), of galaxy 
clustering on large scales (e.g. Cole et alT 20051 ) , and of 
the currently accelerated expansion o f the Universe (e.g 
IClocchiatti et al.ll2006l : IGuy et al.|[201ol ) have confirmed the 
"Lambda cold dark matter " (ACDM) model , first explored 
theoretically 25 years ago l|Davis et al.lll985l ). as the stan- 
dard model of cosmogony. These observations probe a large 
range of scales, from ~ IGpc to ~ lOMpc. On smaller scales, 
where the distribution of dark matter is strongly non lin- 
ear, observational tests of the model are more complicated 
because of the complexity added by galaxy formation pro- 
cesses. However, it is precisely on these scales that the nature 
of the dark matter may be most clearly manifest. For exam- 
ple, if the dark matter is made of warm, rather than cold 
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particles, free streaming in the early universe would have 
erased primordial fluctuations below a scale that depends 
on the mass of the dark matter particle but could be of or- 
der 10^ - 10^°Mq. These mass scales correspond to dwarf 
galaxies and so, in principle, the abundance and properties 
of dwarf galaxies could encode information about the nature 
of the dark matter. 

The validity of the ACDM model on galactic and 
subgalactic scales has been a sub ject of debate for many 
years. Initially Kly pin et all (|l999l ') and lMoore et all (Il999h 
pointed out a large discrepancy between the number of dark 
matter substructures, or subhaloes, that survive inside a 
galactic halo and the number of satellites that are observed 
around the Milky Way. This so-called 'satellite problem' is 
often interpreted as indicating that the model requires most 
of the subhaloes to contain no visible satellite. This aspect 
of the problem, however, is readily solved by invoking the 
known physics of galaxy formation, particularly early reion- 
ization of the intergalactic medium and supernovae feed- 
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back, which inevitably inhibit the formation of stars in 
small mass haloes. Detailed models that reconcile theory 
and observations in this way date back to the early 200 0s 
iBullock et al.ll2000l : iBenson et al.|[2002l : ISomervillel(2002l '). 

The paucity of observed bright satellites, however, is 
only one as pect of the satellite prob lem. As already em - 
phasized bv iKlvpin et al.1 (|l999l ) and iMoore et al.1 (l999l ). 
there is a problem not only with the abundance of satel- 
lites, but also with their distribution of circular velocities. 
In a halo like that of the Milky Way, the ACDM model pre- 
dicts the existence of several subhaloes with maximum cir- 
cular velocities, ymajEl! in excess of ~ 40km s^^. Using the 
high-resolut ion simulations of g al actic haloes of the A quar- 
ius project (|Springel et allbOOSl L lStrigari et al.1 (|2010l ) have 
recently demonstrated that it is possible to find ACDM 
subhaloes that accurately match the observed stellar kine- 
matics of the five well-studied satellites of the Milky Way. 
The best fits, however, invariably have VInax;$40kms~^. [The 
Strigari et al. sample excludes the Large and Small Magel- 
lanic Clouds (LMC and SMC) which reside in more massive 
haloes, and Sagittarius which is currently being disrupted.] 

The discrepancy between the predicted and inferred 
dis tributions of Vmax val u es ha s recently been highlighted 
by iBovlan-Kolchin et al.l (|201ll ). Using also the Aquarius 
haloes, as well as the Via Lactea simulations (|Madau et al.l 
I2OO 8I). they show explicitly that the simulated haloes typ- 
ically contain a few subhaloes which are too massive and 
too dense (as indicated by their value of Vmax/''max) to host 
any of the observed satellites. If such objects existed in the 
Milky Way, they would have to be empty of stars despite 
their mass. This seems very unlikely so, unless the Milky 
Way is atypical, there is an apparent discrepancy between 
model and observations. 

That the Milky Way is not typical of isolated galaxies of 
similar luminosit y and colour has recently been established 
from SDSS data. IlTu et al.l 120111 ) have shown that only 3.5 
per cent of such gala^xies have 2 sa t ellites as bright as the 
Magellanic Clouds, while ICuo et al.l (120111 1 have shown that 
the luminosity function of the bright [My < —14) Milky 
Way satellites has about tw ice the amplitude of the mean 
for similar galaxies (see also lLares et aLlboill ). While these 
measurements show that the Milky Way is not an average 
galaxy, it is not at present possible to compare the distri- 
bution of Vmax of its Satellites with that of similar galaxies 
directly. However, an indirect probe of this distribution can 
be constructed by combining N-body s imulations with a s ub- 
halo abundance matching procedure (|Busha et al.ll201ll ). 

In this paper we explore whether an alternative hypoth- 
esis for the nature of the dark matter, a warm rather than 
a cold particle, can provide a better match to the inferred 
distribution of satellite circular velocities or masses. Specif- 
ically, we test a model in which the dark matter is one 
of the particles predict ed by the 'neutrino minimal stan- 
dard model (i/MSM)' o f lAsaka fc Sh aposhnikov (2005') and 
iBovarskv etal] (|2009al '). In this model there is a triplet of 
sterile neutrinos, the lightest of which could become non- 



The circular velocity is given by y = (GA/(< r)/ry^^, where 
M is the mass enclosed within radius r and G is the universal 
gravitational constant; the value of r at which the maximum of 
this function, Vmax, occurs is denoted by rmax 



Name 


nip [Mq] 


^•200 [kpc] 


M200 [Mq] 


Ns 


Aq-A2 


1.370 X 10* 


245.88 


1.842 X 10^2 


30177 


Aq-A3 


4.911 X 10* 


245.64 


1.836 X 10^2 


9489 


Aq-AW2 


1.370 X 10* 


242.87 


1.775 X 10^^ 


689 


Aq-AW3 


4.911 X 10* 


242.98 


1.778 X 10^2 


338 


Aq-AW4 


3.929 X 10^ 


242.90 


1.776 X 10" 


126 



Table 1. Basic parameters of the simulations analysed in this pa- 
per. The top two simulations are taken fro m the Aquarius sampl e 
of CDM dark matter haloes published in ISpringel et afl l[200^. 
The simulations are of a single halo, Aq-A, at different numeri- 
cal resolutions. The bottom three are WDM counterparts to the 
CDM simulations, as described in the main text. The second to 
fifth columns give the particle mass (nip), the radius of the sphere 
of density 200 times the critical density (r2oo)j the halo mass 
within r200 {M200) and the number of subhaloes within the main 
halo (Ns). The smallest subhaloes, determined by SUBFIND, con- 
tain 20 particles. 



relativistic at a redshift of ~ 10 , have a mass of ^ 2keV, 
and behave as warm dark matter (WDM). This model is 
consistent with astrophysical and particle physics data, in- 
cludin g constraints on neutri no masses from the Lyman-a 
forest jBovarskv eraLll2009bl ). 

To investigate this WDM model we have resimulated 
one of the Aquarius N-body haloes (Aq-A) with the power 
spectrum suppressed at small scales, as expected in the 
WDM case. N-body simulations of galactic and cluster 
WDM haloes were first carried out i n the early 2000s 
llColm. Avila-Reese. fc Valenzuelal |200GI : iBode et al.l I2OOII : 
iKnebe et al.l |2002| ). These studies found that fewer sub- 
haloes form than in the CDM case and that these tend to 
be less concentrated than their CDM counterparts. Quali- 
tatively, we find similar results but the conclusions of these 
early simulations are difficult to interpret because, as we 
shall see later, the sharp cutoff in the power spectrum gives 
rise to the formation of a large nu mber of artificial haloe s 
that are purel y numerical in origin (I Wang fc White! l2007l ). 
More recently, iMaccio fc FontanotI (|2010l ) carried out new 
simulations of WDM models and found that the luminosity 
function of satellites can be reproduced in these models just 
as well as it can in the CDM case. 

Our simulations have orders of magnitude higher res- 
olution than previous ones, enough to investigate reliably 
the inner structure of the galactic subhaloes that are po- 
tential hosts of the dwarf satellites. Furthermore, we carry 
out convergence tests of our results and develop a method 
for distinguishing genuine WDM haloes from the spurious 
objects that inevitably form in simulations of this kind. We 
describe our simulations in Section [5] present our results in 
Section |3l and conclude in Section |4l 



2 THE SIMULATIONS 

To compare the properties of subhaloes in Milky Way mass 
haloes in CDM and WDM universes, we have assembled a 
sample of five high resolution simulations of galactic mass 
haloes. All the simulations have the same basic cosmological 
parameters: in units of the critical density, a total matter 
density, Q,rn = 0.25 and a cosmological constant, ~ 0.75. 
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The linear power spectrum has a spectral index Us = 1 and is 
normalised to Riv e ag = 0.9, with Hp = 100/ikms"^Mpc"^ = 
73kms-^Mpc-^ (jSoringel et al.ll2008 !). Fl 

We have take n two simulations fr om the Aquarius 
project described in lSpringel et al.l (|2008h . both of the same 
halo, Aq-A, but of different resol ution, corresponding to lev- 
els 2 and 3 in the notation of ISpringel et all (l2008h . The 
higher resolution, level 2, simulation has more than a hun- 
dred million particles within r2oo, the radius of a sphere 
about the halo centre, encompassing a mean density of 200 
times the critical density. The level 3 simulation has 3.6 
times fewer particles. In both cases, the mass of the halo 
within r2oo is about 1.8 x lO^^M©, whic h is consistent with 
the estimated i nass of the Milky Way jU fc Whitd l2008l : 
IXue et al.ll200i : lGnedin et al.ll2O10D . The basic properties of 
these haloes are given at the top of Table [T] Substructures 
were identified using the SUBFIND algorithm ( Springcl ct al] 
[2003) to find gravitationally bound subhaloes within them. 

We created three WDM counterparts to the CDM 
haloes by running new sim ulations using the s ame code 
and numerical parameters as 'Sprin gel et al.l (|2008l ) but with 
WDM initial conditions. The WDM initial conditions were 
created keeping the same phases and the same unperturbed 
particle positions as in the CDM case, but using a WDM 
matter power spectrum instead to scale the amplitudes of 
the fluctuations. The linear matter power spectrum for both 
the CDM and WDM simulations is shown in Fig.[l]with solid 
lines adopting an arbitrary normalisation at large scales. 

The WDM power spectrum has a strong cut off at high 
wavenumbers due to the free streaming of the WDM par- 
ticles. In an unperturbed universe at the present day the 
typical velocities of WDM particles are only a few tens of 
ms~^. This implies that the particles ceased to be relativis- 
tic after a redshift of z ~ 10^, well before the end of the 
radiation-dominated era, as suggested by the word 'warm'. 
Fig. [2] illustrates the free streaming of a typical WDM par- 
ticle over cosmic time. The area under the curve is the 
comoving distance traveled. It is evident that the WDM 
particle travels the greatest comoving distance during the 
radiation-dominat ed era after it has become nonrelativistic 
l|Bode et al.ll200ll '). Over the duration of the N-body simu- 
lation, which starts at z = 127, a particle typically travels a 
distance of around 14 kpc, which is small compared to the 
total distance from early times of 400 kpc. For comparison, 
the mean interparticle separation for the high resolution re- 
gion in our highest resolution simulation is 7.4 kpc, similar 
to the free-streaming distance traveled by the particles after 
2 — 127. This means that the effects of streaming during the 
simulation are small, and only affect scales that are barely 
resolved in our simulations. For this reason we chose to set 
the particle velocities in the same way as in the CDM case, 
where the particle velocity is a function of the unperturbed 

^ Although this set of parameters is discrepant at about the 
3(7 level with the latest const raints from microwav e background 
and large-scale structure data jKomatsu et al.ll201ll) . particularly 
with the values of erg and Us, the differences are not imp ortant 
for our purposes. For example, iBovlan-Kolchin et al.l l l201ll'l show 
that the structure of Aquarius subhaloes is statistically similar to 
that of subhaloes in the Via Lactea simulations whic h assume a 
value of (Tg = 0.74, lower than that of lKomatsu et al] l(2011i h and 
a spectral index of 0.95. 




Iog,o(k/h Mpc-') 

Figure 1 . The solid lines show the linear power spectra (from 
CMBFAST; ISeliak fc Zaldarriagal Il996h used for the two simula- 
tions. Black is the original, CDM Aq-A spectrum, and red is 
that of Aq-AW. The vertical dashed line marks the peak of the 
WDM spectrum peak. The arrow marks the Nyquist frequency of 
the level 2 simul ations. The dashed red curve corresponds to the 
M2L25 model of llBovarskv et al .I2OO91J) which is almost identical 
to the solid red curve for scales below fc ^ 10 ?tMpc~^. 

comoving position of a particle and is determined solely by 
the matter fluctuations. 

The WDM matter power spectrum we assume has a 
shape characteristic of a 'thermal relic' (Bode et al. 2001). 
However our WDM matter power spectrum is also an ex- 
cellent fit for scales below fe 10 /iMpc ~^, to the mat- 
ter pow er spectrum of the M2L25 model of iBovarskv et al] 
l|2009bl ). which is shown as a dashed line in Fig. [1] At 
fc = 10 h'mathrmMpc~^ the power in both WDM curves 
is a factor three below that of CDM and falls away very 
rapidly beyond here in both models. The M2L25 model cor- 
responds to a resonantly produced 2keV sterile neutrino 
with a highly non-equilibrium spectrum of primordial ve- 
locities. The mo del is only just consist ent with astrophys- 
ical constraints ()Bovarskv et al.l l2009bl i and so maximizes 
the differences between the substructures in the CDM and 
WDM haloes, both in their internal structure and in their 
abundance. 

For wavenumbers below the peak at 4.5/iMpc~^ the lin- 
ear WDM power spectrum is well approximated by the prod- 
uct of the linear CDM power spectrum times the square of 
the Fourier transform of a spherical top-hat filter of unit 
amplitude and radius 320 kpc, or equivalently, containing a 
mass of 5 X IO^Mq at the mean density. 

Images of the CDM and WDM haloes are shown in 
Fig. |3l As shown in Table [TJ the mass of the main halo in 
the WDM simulation is very similar to that of the CDM 
halo, just a few per cent lighter. However, the number of 
substructures in the WDM case is much lower, refiecting the 
fact that the small scale power in these simulations is greatly 
reduced. Some of the largest subhaloes can be matched by 
eye in the images of the two simulations. 

ISpringel et al.l l|2008h showed that it is possible to make 
precise matches between substructures at different resolu- 
tions for the Aq-A halo, allowing the numerical convergence 
of properties of substructures to be checked for individual 
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Logi„(l+z) 

Figure 2. The free streaming comoving distance traveled per log 
interval o{ 1 + z, where z is redshift, for a WDM particle with 
a fiducial velocity of 24 m s~^ at the present day. The dashed 
vertical line marks the redshift of matter-radiation equality. The 
dotted vertical line indicates the start redshift of the WDM sim- 
ulations. 



substructures. For this paper, we have found matches be- 
tween subhaloes in the Aq-AW2, Aq-AW3, and Aq-AW4 
simulations. We make these matches at the epoch when 
the subhaloes first have a mass which is more than half the 
mass they have at the time when they first infall into the 
main halo (which is very close to the maximum mass they 
ever attain). At this epoch it is relatively easy to match the 
largest substructures in these three simulations as the corre- 
sponding objects have very similar positions, velocities and 
masses. 

The number of subhaloes that can be matched between 
the two WDM simulations is much smaller than that be- 
tween the corresponding CDM simulations, and is also a 
much smaller fraction of the total number of subhaloes iden- 
tified by SUBFIND. The majority of substructures identified 
in the WDM simulations form through fragmentation of the 
sharply delineated filaments characteristic of WDM simu- 
lations and do not have counterparts in the simulations of 
different resolution. The same phenomenon is seen in hot 
dark matter simulations and is numerical in origin, occur- 
ring along t he filaments on a scale matching the interparticle 
separation (|Wang fc Whitell2007l ). This artificial fragmenta- 
tion is apparent in Fig. (3] 

We will present a detailed description of subhalo match- 
ing in a subsequent paper but, in essence, we have found 
that matching subhaloes works best when comparing the 
Lagrangian regions of the initial conditions from which the 
subhaloes form, rather than the subhaloes themselves. We 
use a sample of the particles present in a subhalo at the 
epoch when it had half of the mass at infall to define the 
Lagrangian region from which it formed. We have devised 
a quantitative measure of how well the Lagrangian regions 
of the substructures overlap between the simulations of dif- 



ferent resolution, and select as genuine only those subhaloes 
with strong matches between all three resolutions. We find 
that these criteria identify a sample of fifteen relatively mas- 
sive subhaloes with mass at infall greater than 2 x 10^ M©, 
together with a few more subhaloes with infall mass below 
lO^M©. This sample of fifteen subhaloes includes all of the 
subhaloes with infall masses greater lO^M©. 

We have also found that the shapes of the Lagrangian 
regions of spurious haloes in our WDM simulations are typ- 
ically very aspherical. We have therefore devised a second 
measure based on sphericity as an independent way to reject 
spurious haloes. All fifteen of the massive subhaloes identi- 
fied by the first criterion pass our shape test, but all but one 
subhalo with an infall mass below lO^'Mo are excluded. For 
the purposes of this paper we need only the 12 most massive 
subhaloes at infall to make comparisons with the Milky Way 
satellites. 

For both our WDM and CDM catalogues, we select a 
sample made up of the 12 most massive subhaloes at in- 
fall found today within 300 kpc of the main halo centre. In 
the Aq-AW2 simulation these subhaloes are resolved with 
between about 2 and 0.23 million particles at their maxi- 
mum mass. We use the particle nearest the centre of the 
gravitational potential to define the centre of each subhalo 
and hence determine the values of Knax and rmax defined in 
Section 1. 



3 RESULTS 

In this section, we study the central masses of the sub- 
structures found within 300 kpc of the centres of the CDM 
and WDM Milky Way like haloes. These results are com- 
pa red with the masses within the half-light r a dii, in ferred 
by IWalker et all (|2009l . |2010| ) and IWolf et al.1 (|2010D from 
kinematic measurements, for the 9 bright [Lv > IO^Lq) 

Milky Way dwarf spheroid al galaxies. 

Following the study bv lBovlan-Kolchin et al.l (|201ll ). in 
Fig. U we plot the correlation between Vmax and Tmax for 
the subhaloes in Aq-AW2 and Aq-A2 that lie within 300kpc 
of the centre of the main halo. Only those WDM subhaloes 
selected using our matching scheme are included, whereas 
all Aq-A2 subhaloes are shown. Th e CDM subhaloes are 
a sub set of those shown in fig. 2 of iBovl an-Kolchin et al.l 
(|201lD . and show Vmax values that are typically ~ 50 per 
cent larger than those of WDM haloes with a similar r^ax. 
By assuming that the mass density in the subhaloes contain- 
ing the observed dwar f sphe r oidals follows an NFW profile 
iNavarro et al.lll996bl . Il997l ). iBovlan-Kolchin et al.l (|201lh 
found the locus of possible (rmax, Knax) pairs that are con- 
sistent with the observed half-light radii and their enclosed 
ma sses. This is represented by the shaded region in Fig. [l] 
As iBovlan-Ko lchin et al. ( 2011) observed with their larger 
sample, several of the largest CDM subhaloes have higher 
maximum circular velocities than appears to be the case for 
the Milky Way bright dwarf spheroidals. By contrast, the 
largest WDM subhaloes are consistent with the Milky Way 
data. 

Rather than assuming a functional form for the mass 
density profile in the observed subhaloes, a more direct ap- 
proach is to compare the observed masses within the half- 
light radii of the dwarf spheroidals with the masses within 
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Figure 3. Images of the CDM (left) and WDM (right) level 2 haloes at 2: = 0. Intensity indicates the line-of-sight projected square 
of the density, and hue the projected density-weighted velocity dispersion, ranging from blue (low velocity dispersion) to yellow (high 
velocity dispersion). Each box is 1.5 Mpc on a side. Note the sharp caustics visible at large radii in the WDM image, several of which 
are also present, although less well defined, in the CDM case. 
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Figure 4. The correlation between subhalo maximum circular 
velocity and the radius at which this maximum occurs. Sub- 
haloes lying within 300kpc of the main halo centre are in- 
cluded. The 12 CDM and WDM subhaloes with the most mas- 
sive progenitors are shown as blue and red filled circles respec- 
tively; the remaining subhaloes are shown as empty circles. The 
shaded area represents the 2a confidence region for possible hosts 
of the 9 bright Milky Wa y dwarf spheroidals determined by 



ot tne V brignt IViUky Wa y 
iBovlan-Kolchin et al.l boill) . 



the same radii in the simulated subhaloes. To provide a fair 
comparison we must choose the simulated subhaloes that 
are most likely to correspond to those that host the 9 bright 
dwarf spheroidals in the Milky Way. As stripping of sub- 
haloes preferentially removes dark matter relative to the 
more centrally concentrated stellar component, we choose to 



associate final satellite luminosity with the maximum pro- 
genitor mass for each surviving subhalo. This is essentially 
the mass of the object as it falls into the main halo. The 
smallest subhalo in each of our samples has an infall mass 
of 3.2 X lO^M© in the WDM case, and 6.0 x IO^Mq in the 
CDM case. 

The LMC, SMC and the Sagittarius dwarf are all 
more luminous than the 9 d warf spheroidals considered by 
iBovlan-Kolchin et al.l (1201 ll ) and by us. As noted above, the 
Milky Way is exceptional in hosting galaxies as bright as 
the Magellanic Clouds, while Sagittarius is in the process of 
being disrupted so its current mass is difficult to estimate. 
Boylan-Kolchin et al. hypothesize that these three galaxies 
all have values of Vmax > 60kms~^ at infall and exclude sim- 
ulated subhaloes that have these values at infall as well as 
Knax > 40kms~^ at the present day from their analysis. In 
what follows, we retain all subhaloes but, where appropri- 
ate, we highlight those that might host large satellites akin 
to the Magellanic Clouds and Sagittarius. 

The circular velocity curves at z = for the 12 sub- 
haloes which had the most massive progenitors at infall are 
shown in Fig. [S] for both WDM and CDM. The circular 
velocitie s within the half-l ight radius of the 9 satellites mea- 
sured by I Wolf et al.l (|2010l ) are also plotted as symbols. Leo- 
II has the smallest half-light radius, 200pc. To compare 
the satellite data with the simulations we must first check 
the convergence of the simulated subhalo masses within at 
least this radius. We find that the median of the ratio of the 
mass within 200pc in the Aq-W2 and Aq-W3 simulations is 
W2/W3 - 1.22, i.e., the mass within 200pc in the Aq-W2 
simulation has converged to better than ~ 22%. 

As can be inferred from Fig. O the WDM subhaloes 
have similar central masses to the observed satellite galax- 
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Figure 5. Circular velocity curves for the 12 CDM (left) and 
WDM (right) subhaloes that had the most massive progenitors. 
The 3 red curves represent subhaloes with the most massive pro- 
genitors, which could correspond to those currently hosting coun- 
terparts of the LMC, SMC and the Sagittarius dwarf. The 9 black 
curves might more fairly be compared with the data for the 9 
bright dwa r f sphe roidal galaxies of the Milky Way considered by 
IWolf et al.l 1 I2OICII ). Deprojected half-light rad ii and their corre- 
sponding half-light masses, as determined bv IWolf et al. 1 1I2OIOI ') 
from line-of-sight velocity measurements, are used to derive the 
half-light circular velocities of each dwarf spheroidal. These veloc- 
ities and radii are shown as coloured points. The legend indicates 
the colour coding of the different galaxies. 



ies, while the CDM subhaloes are almost all too massive at 
the corresponding radii. The CDM subhaloes have central 
masses that are typically 2-3 times larger than the Milky 
Way satellites. There is one CDM subhalo that lies at lower 
masses than all 9 dwarf spheroidals, but this had one of the 
three most massive progenitors and has been almost com- 
pletely destroyed by tidal forces. 

Fig. |4] and [5] show that the WDM subhaloes are less 
centrally concentrated than those in the corresponding CDM 
halo. Concent rations typically reflect the epoch at which th e 
halo formed (jNavarro et all Il996bl . 1 19971 : lEke et al.l l200lh . 
To investigate systematic differences in the formation epoch 
of the WDM and CDM subhaloes in our sample, we must 
choose a suitable definition of formation time. Since we are 
considering only the central mass, and we do not wish to 
introduce scatter in any correlation by using subhaloes that 
may have been stripped, we define the formation time as 
the first time at which the total progenitor mass exceeds the 
mass within 1 kpc at infall. The correlation of this redshift 
with the mass within 1 kpc at infall is shown in Fig.[6]for the 
12 most massive WDM and CDM progenitors that survive to 
z = as distinct subhaloes. Evidently, the proto subhaloes 
that form later, which are generally WDM not CDM ones, 
have the lowest central masses. The mean difference between 
the top 12 WDM and CDM proto-subhalo masses within 1 
kpc is approximately a factor 2. 

Because of their later formation time, the infalling 
WDM subhaloes already have lower central masses than 
those falling into the corresponding CDM haloes. As their 
mass is less centrally concentrated, the WDM subhaloes are 
more susceptible to stripping. While this is most impor- 
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Figure 6. The correlation between subhalo central mass at infall 
and the redshift of formation, 2formi defined as the redshift at 
which the total mass of each proto subhalo first exceeded this 
value. Central mass is defined within 1 kpc, and CDM and WDM 
results arc shown with blue and red symbols respectively. 



tant in the outer regions of the subhaloes, the mass profiles 
in Fig. [5] show that the inner regions of some of the sub- 
haloes have also endured significant depletion since infall. 
Fig. [7] shows, for both WDM and CDM subhaloes, the ra- 
tio, Af2=o(< '")/Afinfaii, of the present day mass contained 
within r = 0.5, 1 and 2 kpc to the mass at infall, as a 
function of the central mass at infall at the chosen radius. 
On average, the median mass at infall for WDM is lower 
by ^ 0.15 dex than the corresponding mass for CDM. One 
subhalo gains mass between infall and z = Q because it ac- 
cretes another subhalo. While there is a large scatter among 
the different subhaloes, with some having lost the majority 
of their central mass since infall, no significant systematic 
difference between WDM and CDM subhaloes is apparent. 
This implies that the reason why the WDM subhaloes pro- 
vide a better fit to the half-li ght masses of the 9 Milky Way 
dwarf spheroidals studied bv iWolf et al.l (|2010l ) is not excess 
stripping but the later formation time, and correspondingly 
typical lower concentration, of the WDM proto subhaloes 
compared to their CDM counterparts. 



4 DISCUSSION AND CONCLUSIONS 

The properties of the satellite galaxies of the Milky Way 
have posed a longstanding puzzle for CDM theories of galaxy 
formation. Two aspects of this puzzle have reportedly been 
separately and independently solved. One is the luminos- 
ity function of the satellites. The basic idea - the suppres- 
sion of galaxy formation in small haloes by a combination 
of feedback effects produced by the reionization of gas at 
high redshift and supernova hea t ing - was suggested by 



iKauffmann. White, fc Guiderdoni (Il993) and explored thor- 
ough l y in the early 2000s (jBuUock et al.ll2000l : iBenson et al.1 
I2OO2I : ISom ei^iU c 2002) and has been revisited many times 
since then (see iFont et al.ir201ll . and references therein for 
the most recent discussion). The other aspect concern s the 
dynamical state of the satellites. IStrigari et al.l l|2010l ) have 
shown that there exist subhaloes in the Aquarius CDM sim- 
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Figure 7. The variation with subhalo mass at infall of the ratio 
of the present day mass to the infall mass contained within 500pc, 
Ikpc and 2kpc. Data are shown for the 12 subhaloes identified at 
z = which had the most massive progenitors, with CDM in blue 
and WDM in red. The symbol type denotes the radius interior 
to which the central mass is being measured and large symbols 
show the medians of the corresponding distributions. We find no 
systematic differences between the CDM and WDM subhalo mass 
ratios. 



ulations that fit the steUar spectroscopic data for the weU- 
studied satelhtes extremely well. 

There is a third aspect to the puzzle, however, that 
has not yet been fully addressed and this is whether the 
CDM models that account for the satellite luminosity func- 
tion also account for the satellites' internal dynamics. In 
other words, do the models assign the correct luminosities 
to subhaloes with the correct dynamics? At face value, the 
answer seems to be 'no'. Th is is already evident in the anal- 
ysis of IStrigari et al.l (|2010l ) in which the best fit dynami- 
cal models imply velocity dispersions (or equivalently Vmax 
values) for the brightest dwarf spheroidals that are smaller 
than the velocity dispersions of the largest subhaloes. It 
is this discrepancy that ha s recently been highlighted by 
iBovlan-Kolchin et al.l (120111 ). 

In this paper, we have compared a high resolution N- 
body simulation of one of the Aquarius galactic haloes with 
a WDM counterpart. The initial conditions for both had 
the same phases and the same unperturbed particle posi- 
tions. For the WDM simulation we chose a form of the power 
spectrum corresponding to on e of the models discussed b y 
lAsaka fc Shaposhnikov (|2005l ) and iBovarskv et al.l (|2009al ). 
in which the dark matter is a sterile neutrino with mass ~ 
2keV, just consistent w ith various astrophysical constraints 
l|Bovarskv et al.ll2009bl ). The suppression of the power spec- 
trum at masses below ~ IO^^Mq delays the formation of the 
haloes that will end up hosting the satellites and, as we have 
shown, this lowers their concentration compared to that of 
the corresponding CDM haloes. This is enough to reconcile 
the dynamics of the subhaloes with the data. 

While a WDM model naturally produces haloes that 
are less concentrated than their CDM counterparts, this is 
only one possible solution to the puzzle. Other forms of 
dark matter such as 'meta-CDM' resulting from the de- 
cay of cold thermal relics could produce a similar outcome 



l|Strigari et al.ll2007l ). Also, it must be borne in mind that 
the values of Knax for Milky Way satellites are not di- 
rectly measured but inferred by making assumptions about 
their dynamical state. If some of these assumptions are un- 
realistic, this coul d lead to an under estimate of the val- 
ues of Knax (e.g. IStoehr et al.l |2002| ). Another possibil- 
ity is that the satellite population of the Milky Way is 
not typical of the average to which the model pre d iction s 
apply. It h as re cently been show n by Liu et ahl l|201ll ) , 
ICuo et al.1 (|201ll ') and iLares et al] (|201ll ) that the bright 
end of the Milky Way satellite luminosity function is differ- 
ent from the average. Finally, we cannot exclude the pos- 
sibility that baryonic processes occurring during the for- 
mation of satellite galaxies in the CDM cosmogony might 
have lowered the concentr ation of haloes , for ex ample, in 
the manner suggested by 'Navarro et al.' (l996a!). Recent 
simulations iRcad fc G ilmorc 2005; Mashchcnko ct al. 20081; 
iGovernato et al.ll20ld ) suggest that these processes could be 
important although it remains to be seen if they are enough 
to reconcile the CDM model with the dynamics of the Milky 
Way satellites. 
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